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We study the quantum theory of a simple general relativistic quantum model of two coupled 
harmonic oscillators and compute the two-point function following a proposal first introduced in the 
context of loop quantum gravity. 



I. INTRODUCTION 

The quantum dynamics of a general relativistic system is not governed by a genuine Hamiltonian. Instead 
the evolution is fixed by the first class constraints of the theory, among which is the Hamiltonian constraint, 
and therefore evolution does not refer to an external (classical) time parameter - as is the case in standard 
quantum mechanics - but corresponds to the unfolding of gauge transformations. The observables of the 
theory are gauge invariant quantities represented by operators commuting with the first class constraints 
and in particular with the Hamiltonian constraint and therefore are all constants of motion. This is at the 
heart of the so called "problem of time" for general relativistic systems: How can we reconcile a timeless 
quantum evolution with the experimental evidence of the flow of time? A related question is how can we 
extract from the quantum system a classical limit where the dynamical evolution is described with respect 
to a parameter playing the role of time. 

In order to study such issues we examine a simple model of two harmonic oscillators characterized by an 
hamiltonian constraint. Such a model was first considered by Rovelli in [lj] and its quantum dynamics was 
examined in details in [2]. In particular in this paper we compute and study the two-point function of this 
model, following a proposal introduced in [3] to compute the graviton propagator in loop quantum gravity. 
In Sec. II we describe the classical and quantum dynamics of the system of two harmonic oscillators with 
an Hamiltonian constraint, and we give the expression of the projector operator from the kinematical to the 
physical Hilbert space of the quantum theory. Coherent states are introduced in Sec. III. They will play a 
key role in the following. Finally we compute explicitly in Sec. IV the two-point function of the model and 
comment some of its properties. 

II. DOUBLE HARMONIC OSCILLATORS 

Let us consider a classical system of two harmonic oscillators (of the same frequency and the same mass 
equal to 1) with positions q a , q b and momenta p a ,Pb- The dynamics is governed by the hamiltonian constraint 

H(q a) q b , Pa)Pb ) = -(q 2 a + ql+pl+p 2 b )-M ~ 0, (1) 

where M is a real positive constant characterizing the system. The extended configuration space is C = M 2 
coordinatized by the observables quantities q a and q b , called partial observables 1 . The function H is defined 
on the cotangent space T*C = M 4 . 
The action of the system is 

S = f dr ( Pa q a + Pb q b - iV~ {q 2 a + q 2 b + p\ + p\ - 2M) J , (2) 
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where a dot denotes the derivative with respect to t and N is a lagrangian multiplier. The variation of 
the action with respect to the variables q a (T~),p a (T),qb(T),pb(T) and N(t) gives the canonical equations of 
motion and the constraint equation {1} . The action is invariant under reparametrization of the parameter r 
labeling the motion. 

The solution to the equations of motion can be parametrized by this non-physical parameter r 

q a (T) = A sm(r + a ), (3) 
q b (r) = B sin(r + <&,). (4) 

This motion is an ellipse of radii A and B and inclination A</> = <f> a — <f>b in the extended configuration space 

C, the plane (q a , with (A, B) £ 0, y/2M and (</>, 0') € [0, 27r]. The values of the constants A, £? are fixed 

by the constraint {J) to satisfy A 2 + B 2 — 2M. Due to the reparametrization invariance, the parameter r 
is not connected with observability and the physical predictions of the system regard relation between the 
partial observables a and b, not their dependence on t. A gauge transformation changes the phases 4> a and 
4>b but not A(f>. Therefore the physical phase space is coordinatized by the two quantities A and A(j> which 
are constant along the motions, and identify the ellipse. 

The canonical quantization of the system is straightforward: The partial observables q a , qb and p a ,Pb are 
promote to self-adjoint operators acting on the kinematical Hilbert space /C = dq a dqb] as multiplicative 

operators and derivative operators respectively, and satisfying the usual commutation relations. The physical 
Hilbert space TL of the theory is the space of solutions of the equation defining the dynamics, i.e. the Wheeler- 
DeWitt equation: 

d 2 d 2 \ 
-Qq2 -gj2 +92+3?-2Mj^(g a ,<ft) = 0, (5) 

where we have take K=l. We introduce the creation and annihilation operators for the two oscillators, with 
the usual commutators [a, a'] = [b,w] — 1, [a,b] = [a,b'] = 0, and the energy eigenstates \n a ,n b ) where n a 
and rib are the energy levels of the oscillators q a and qb respectively. The physical Hilbert space is spanned 
by the vectors satisfying the equation M = n a + n b + 1. It is convenient to introduce the quantum number 
to = h{n a — nb) that runs from m = — j to to = j, with j = (M — l)/2. So that TL is spanned by the (2j + 1) 
states \m) = \j — m,j + in). The representation of the states in the space of coordinates is 

i> m {a,b) = (2 2j '7r(i + m)!(i-m)!)" 1/2 ff J+m (a) J H i _ m (6)e-^ i , (6) 

where H n (q) id the n-th Hermite polynomial. The dynamics can be obtained from an orthogonal projection 
operator P from the kinematical Hilbert space /C into the physical Hilbert space H, P : JC — > TL, defined by 

P = J dre~" H = J2 HH, (7) 

m=-j 

and the propagator K of the model is the integral kernel of P. The properties and the classical limit of the 
propagator K have been studied in 

III. COHERENT STATES 

The coherent states in the kinematical Hilbert space K can be written in terms of the energy eigenstates 

\a^) = e-^-^^^^\n a ,n b ), (8) 

where the parameters a, (3 are functions of the classical positions (q a , qb) and the classical momenta {p a ,pb): 
a = q a + ip a and (3 — qb + ipb- To obtain the coherent states in the physical Hilbert space we use the 
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projector P: 

\a,/3)phys = P\ot,0) 

_ -\ a \ 2 /2-\P\ 2 /2 



A 2 /2-\f}\ 2 /2 ST 



/ --\n 

*=-,■ v (i - m ) ! (i + m ) ! 



(9) 



The wave function of this coherent state is peaked around the classical trajectory determined by the initial 
conditions {q a ,qb,Pa,Pb), namely around the ellipse with parameter A and A<fi determined by a and (3. The 
normalized coherent state is 

\a,P)ph v 7 = , ) = =; \a,P) P hys 



Following the interpretation of general relativistic quantum mechanics given in the scalar product between 
two coherent states represents a physical transition amplitude, 

(a'a + 0'3) 2j 

ph„ a \ ,H, ph y S (|o'|2 + |/3'|2) J (|o| 2 + |/3| 2 ) J 

This product is in general different from zero, revealing the existence of some overlap between coherent 
states peaked around different classical solutions. However the modulus square of expression (fTTj) . having 
the interpretation of the transition probability form the state to the state \a' , /?')ph™> results to 

be maximal when a = a' and j3 = /?'. In Q Ashworth expresses the scalar product (fTT|l in terms of the 
quantity £ = a//3, and shows that this scalar product, seen as a function of, say, £' for £ fixed, is a Gaussian 
centered on £ with width equal to (2j)~ 1 . Therefore in the classical limit where 2j is large the transition 
between coherent states peaked around different classical trajectories is highly suppressed. 



IV. TWO-POINT FUNCTION 



The two-point function for a single (non-relativistic) harmonic oscillator is given by 

G(h,t 2 ) = (0|z(ti):r(i2)|0) = (0\x e -^-*=) x \0), (12) 

where H is the hamiltonian and K = 1. The standard interpretation is that G(t%, t 2 ) represents the probability 
amplitude that a particle created at the instant t\ is measured at the instant t 2 . This definition of the two- 
point function appears to be hardly implementable for a general relativistic quantum mechanical system. The 
difficulties come from the absence of a vacuum state and a time parameter with respect to which evolution 
is defined for such system. How can we define, in a general relativistic quantum mechanical system, the 
distance between the two points appearing in lfl2|) ? How can we give an interpretation in terms of a particle 
propagating between the two points? 

Recently a technique to compute the two-point function in the context of quantum gravity has been 
proposed in 0| . We want to apply this proposal to our system of two harmonic oscillators that share with 
quantum gravity the basic property of not having any background time parameter and no state minimizing 
the energy, i.e. no vacuum state. The ingredients of the proposal are a boundary semiclassical state, the 
propagator of the theory, and operators acting on the boundary state. Following 0| we take as boundary 
state a state living in the Hilbert space Ji* x H, in coordinate representation 

^{qa,qb,q' a ,q' b ) = ^i{q a , Qb) ip2(q' a , (13) 

where ^i{q a ,q b ) = (q a , qb\a\, and ip 2 (q a ,qb) = (q a , qb\&2, fa)^™- The use of coherent state is the 

key idea in order to define the two points appearing in lfl2|) . better the "temporal distance" between the 
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two points. We choose the two coherent states ipi and tp2 as peaked around the same classical trajectory. 
Moreover the parameter ati , f3\ , , /?2 will identify two particular points on this trajectory. For an ellipse 
given by the classical motion l|3|4p . these parameters can be written as 



ai = Ae- 4(ri+< H (14) 

ft = Be-^ 1 ^, (15) 

a 2 = Ae- i( - T * +M , (16) 

/3 2 = Be- i(T2+ *>. (17) 

Therefore the distance between the two points will be measured by the difference of the phase factors, 
At = ti — r 2 . 

The next ingredient is easly obtained, it's the propagator of the model, namely the matrix elements of the 
projector (0 on the physical Hilbert space, K(q a ,q b ,q' a ,q' b ) = (q a ,q b \P\q' a ,q' b ), 

Finally we need the operators acting on the boundary state. For that we simply take the product of the 
positions operators q a and q b . We notice that these operators are not gauge invariant, because they do not 
commute with the hamiltonian constraint. In particular the product q a q b contains a gauge invariant part, 
given in terms of the creation and annihilation operators by atf + a^b, and a non gauge invariant part, 
ab + a^tf . The role of the propagator will be the one to select the gauge invariant part of q a q b . 

We are now ready to compute the two-point function for our general relativistic quantum mechanical 
system following the proposal of [3|. We indicate the two-point function as G(At), where the argument of 
G is connected with the distance of the two points. The expression is 

G(At) = J dq a dq b dq' a dq b ^ 1 (q a ,qb)q a qbK(q a ,q b ,q' a ,q b )q' a q b Jp2(q'a,qb), ( 18 ) 
or in the Dirac notation 

G(Ar)= ^(^fclqaqbPqaqbla!,^};^?- (19) 
With simple algebra we obtain 

(a 2 -a 1 +^/3 1 ) 2 ^ 2 ~~ 



G(Ar) = 2j(2j - 1) 



(M 2 + |ft| a WM 2 + l/W 
Substituting the parameters a and (3 with their expression l|14H17p 



(P2«i) + (a 2 pi) + 2a 2 ai/? 2 pi + 



2J-1 



(.20) 



G(At) = 2j( ^ /2 1} [A 2 B 2 cos 2 A0 + ^--j e-^ AT . (21) 

In the semiclassical limit where 2j ~ M the two-point function takes the form 

G(At) = (A 2 B 2 cos 2 A0 + M) er lMAr . (22) 
The result deserves various comments. 

• The two-point function is a gauge invariant quantity, and therefore can be seen as real physical pre- 
diction of the theory. This can be easly proved noting that a gauge transformation on the system 
transforms the complex quantities as a — > ae l9 and (3 — > fie 10 . In the coherent states these gauge 
transformations appear as an overall phase in front of the vector, |o<i,Pi)p£™ ~~ * e * 2 "' e l a i"^i)ph™- 
Hence the expression (fl9ll of the two-point function is manifestly gauge invariant. In fact the resulting 
formula lf2lj) (as well as the semiclassical limit |22|) ) depends only on the difference of the phases of 
(|14H17p : a gauge transformation simply adds a constant quantity to these phases, and such quantity 
disappears in the difference. 



All the information about the distance of the two points is contained in the boundary state, namely in 
the coherent states determined by parameters l|14H17p . The presence of the propagator in the expression 
of the two-point function has the unique role to select the gauge invariant part of the product of the 
positions operators, q a qt- In contrast with the two-point function for the simple harmonic oscillator 
(fl2l) where the notion of the two points is coded in the propagator exp —iH(t\ — t 2 ), the action of the 
propagator in lfl8|) only guarantees the gauge invariance of the full expression. 
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• The structure of the result f2Tj) (or the semiclassical limit lj22|) ) has some analogies with the two-point 
function of a simple harmonic oscillator (fl2"ll . An easy computation gives (for convenience we explicitly 
write the mass m and the frequency w of the oscillator) 

G(t u t 2 ) = (Olx^M^lO) = J—e-'i^-^. (23) 

2muj 

Hence the two-point function of a simple harmonic oscillator is composed by two factors, a constant 
term depending only on the characteristics of the system, and a complex exponential depending on the 
energy of the propagating particle and the time interval of the propagation. The two-point function for 
our system has a similar structure: Two factors, the first is a constant term function of the properties 
of the system (the quantity M) and of the particular classical trajectory selected by the boundary state 
(the parameters A and A</>), and the second is a complex exponential of the "time interval" separating 
the two points. 

It is worth noticing that the factor |w appearing in the exponential of (|23| represents the energy of 
all the particles involved in the process: In this case it is the energy of the first excited state of the 
quantum harmonic oscillator. In the very same way, the factor 2j in the exponential of the two-point 
function lj2Tj) represents the total energy of the particle taking part in the process: Such energy is fixed 
by the quantum hamiltonian constraint. 

• Another property of the two-point function lfT8|) is to be proportional to the scalar product of the 
coherent states that define the boundary state. In fact the scalar product of the physical states 

&>$ST and \a a ,fo)$F is given by 

™ / °™(a 2 ,/? 2 |a 1 ,/3 1 )^™ = e'^ A \ (24) 
which is exactly the exponential appearing in (|2"Tj) . 



V. SUMMARY 



We have computed the two-point function for the system of two quantum harmonic oscillators which 
dynamics is fixed by an Hamiltonian constraint. The computation was based on the proposal of 0] for the 
two-point function in the context of loop quantum gravity. In our model a key ingredient has been the use of 
coherent states defined in the physical Hilbert space and peaked around classical solutions. The information 
concerning the two points of the correlation function is entirely contained in the boundary state build on 
the tensor product of two coherent states. The notion of time and in particular the notion of the "temporal 
distance" between the two points emerges from an appropriate choice of the boundary state, and is not 
related with the propagator of the system whose role is only to ensure the gauge invariance of the two-point 
function. 
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